\(\int x^m (1-a x)^n (1+a x)^n \, dx\) [986]

Optimal result

Integrand size = 19, antiderivative size = 36 \[ \int x^m (1-a x)^n (1+a x)^n \, dx=\frac {x^{1+m} \operatorname {Hypergeometric2F1}\left (\frac {1+m}{2},-n,\frac {3+m}{2},a^2 x^2\right )}{1+m} \]

[Out]

x^(1+m)*hypergeom([-n, 1/2+1/2*m],[3/2+1/2*m],a^2*x^2)/(1+m)

Rubi [A] (verified)

Time = 0.01 (sec) , antiderivative size = 36, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.105, Rules used = {126, 371} \[ \int x^m (1-a x)^n (1+a x)^n \, dx=\frac {x^{m+1} \operatorname {Hypergeometric2F1}\left (\frac {m+1}{2},-n,\frac {m+3}{2},a^2 x^2\right )}{m+1} \]

[In]

Int[x^m*(1 - a*x)^n*(1 + a*x)^n,x]

[Out]

(x^(1 + m)*Hypergeometric2F1[(1 + m)/2, -n, (3 + m)/2, a^2*x^2])/(1 + m)

Rule 126

Int[((f_.)*(x_))^(p_.)*((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int[(a*c + b*d*x^2)
^m*(f*x)^p, x] /; FreeQ[{a, b, c, d, f, m, n, p}, x] && EqQ[b*c + a*d, 0] && EqQ[n, m] && GtQ[a, 0] && GtQ[c,
0]

Rule 371

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[a^p*((c*x)^(m + 1)/(c*(m + 1)))*Hyperg
eometric2F1[-p, (m + 1)/n, (m + 1)/n + 1, (-b)*(x^n/a)], x] /; FreeQ[{a, b, c, m, n, p}, x] &&  !IGtQ[p, 0] &&
(ILtQ[p, 0] || GtQ[a, 0])

Rubi steps \begin{align*} \text {integral}& = \int x^m \left (1-a^2 x^2\right )^n \, dx \\ & = \frac {x^{1+m} \, _2F_1\left (\frac {1+m}{2},-n;\frac {3+m}{2};a^2 x^2\right )}{1+m} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.06 (sec) , antiderivative size = 38, normalized size of antiderivative = 1.06 \[ \int x^m (1-a x)^n (1+a x)^n \, dx=\frac {x^{1+m} \operatorname {Hypergeometric2F1}\left (\frac {1+m}{2},-n,1+\frac {1+m}{2},a^2 x^2\right )}{1+m} \]

[In]

Integrate[x^m*(1 - a*x)^n*(1 + a*x)^n,x]

[Out]

(x^(1 + m)*Hypergeometric2F1[(1 + m)/2, -n, 1 + (1 + m)/2, a^2*x^2])/(1 + m)

Maple [F]

\[\int x^{m} \left (-a x +1\right )^{n} \left (a x +1\right )^{n}d x\]

[In]

int(x^m*(-a*x+1)^n*(a*x+1)^n,x)

[Out]

int(x^m*(-a*x+1)^n*(a*x+1)^n,x)

Fricas [F]

\[ \int x^m (1-a x)^n (1+a x)^n \, dx=\int { {\left (a x + 1\right )}^{n} {\left (-a x + 1\right )}^{n} x^{m} \,d x } \]

[In]

integrate(x^m*(-a*x+1)^n*(a*x+1)^n,x, algorithm="fricas")

[Out]

integral((a*x + 1)^n*(-a*x + 1)^n*x^m, x)

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 15.68 (sec) , antiderivative size = 206, normalized size of antiderivative = 5.72 \[ \int x^m (1-a x)^n (1+a x)^n \, dx=\frac {a^{- m - 1} {G_{6, 6}^{5, 3}\left (\begin {matrix} - \frac {m}{2} - \frac {n}{2}, - \frac {m}{2} - \frac {n}{2} + \frac {1}{2}, 1 & \frac {1}{2} - \frac {m}{2}, - \frac {m}{2} - n, - \frac {m}{2} - n + \frac {1}{2} \\- \frac {m}{2} - n - \frac {1}{2}, - \frac {m}{2} - n, - \frac {m}{2} - \frac {n}{2}, - \frac {m}{2} - n + \frac {1}{2}, - \frac {m}{2} - \frac {n}{2} + \frac {1}{2} & 0 \end {matrix} \middle | {\frac {e^{- 2 i \pi }}{a^{2} x^{2}}} \right )} e^{- i \pi m} e^{- i \pi n}}{4 \pi \Gamma \left (- n\right )} - \frac {a^{- m - 1} {G_{6, 6}^{2, 6}\left (\begin {matrix} - \frac {m}{2} - \frac {1}{2}, - \frac {m}{2}, \frac {1}{2} - \frac {m}{2}, - \frac {m}{2} - \frac {n}{2} - \frac {1}{2}, - \frac {m}{2} - \frac {n}{2}, 1 & \\- \frac {m}{2} - \frac {n}{2} - \frac {1}{2}, - \frac {m}{2} - \frac {n}{2} & - \frac {m}{2} - \frac {1}{2}, - \frac {m}{2}, - \frac {m}{2} - n - \frac {1}{2}, 0 \end {matrix} \middle | {\frac {1}{a^{2} x^{2}}} \right )}}{4 \pi \Gamma \left (- n\right )} \]

[In]

integrate(x**m*(-a*x+1)**n*(a*x+1)**n,x)

[Out]

a**(-m - 1)*meijerg(((-m/2 - n/2, -m/2 - n/2 + 1/2, 1), (1/2 - m/2, -m/2 - n, -m/2 - n + 1/2)), ((-m/2 - n - 1
/2, -m/2 - n, -m/2 - n/2, -m/2 - n + 1/2, -m/2 - n/2 + 1/2), (0,)), exp_polar(-2*I*pi)/(a**2*x**2))*exp(-I*pi*
m)*exp(-I*pi*n)/(4*pi*gamma(-n)) - a**(-m - 1)*meijerg(((-m/2 - 1/2, -m/2, 1/2 - m/2, -m/2 - n/2 - 1/2, -m/2 -
n/2, 1), ()), ((-m/2 - n/2 - 1/2, -m/2 - n/2), (-m/2 - 1/2, -m/2, -m/2 - n - 1/2, 0)), 1/(a**2*x**2))/(4*pi*g
amma(-n))

Maxima [F]

\[ \int x^m (1-a x)^n (1+a x)^n \, dx=\int { {\left (a x + 1\right )}^{n} {\left (-a x + 1\right )}^{n} x^{m} \,d x } \]

[In]

integrate(x^m*(-a*x+1)^n*(a*x+1)^n,x, algorithm="maxima")

[Out]

integrate((a*x + 1)^n*(-a*x + 1)^n*x^m, x)

Giac [F]

\[ \int x^m (1-a x)^n (1+a x)^n \, dx=\int { {\left (a x + 1\right )}^{n} {\left (-a x + 1\right )}^{n} x^{m} \,d x } \]

[In]

integrate(x^m*(-a*x+1)^n*(a*x+1)^n,x, algorithm="giac")

[Out]

integrate((a*x + 1)^n*(-a*x + 1)^n*x^m, x)

Mupad [F(-1)]

Timed out. \[ \int x^m (1-a x)^n (1+a x)^n \, dx=\int x^m\,{\left (1-a\,x\right )}^n\,{\left (a\,x+1\right )}^n \,d x \]

[In]

int(x^m*(1 - a*x)^n*(a*x + 1)^n,x)

[Out]

int(x^m*(1 - a*x)^n*(a*x + 1)^n, x)